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Recent measurements in the charm sector are reviewed, concentrating on results which are sensitive to New Physics 
effects. The scope of the presentation includes D° — D° mixing searches, a CPT / Lorentz invariance study, and a range 
of searches for rare and forbidden decays. Results from the BaBar, Belle, CDF, CLEO, and FOCUS collaborations 
are presented, including an important first observation. 



1. Introduction 

Presentations on "New Physics" can produce a feel- 
ing of anti-climax, for surely if there were any New 
Physics signals to report, the news would have leaked 
out. One does not expect to hear anything new. This 
talk does contain at least one first observation, how- 
ever, and I'll try to maintain some suspense by not 
mentioning in advance what it is. 

1.1. What I'm Not Talking About 

Today we are at a particular disadvantage because 
the year's biggest charm news is outside the scope of 
the talk. BaBar's discovery 1 of a narrow resonance 
decaying to D s ir° came as a complete surprise — 
apart from the familiar D*, no mesons decaying to 
this final state were foreseen — and led to a flurry 
of speculation that the new D s j(2317) might be an 
exotic meson. The honors were shared among the B- 
factories in an amusing way: CLEO announced the 
discovery 2 of a second state, the D s j (2457), decay- 
ing to -D*7T°; and Belle, as well as confirming these 
results, 3 made the first observation of both states in 
B meson decays. 4 The decay modes and widths of 
the new D s j are consistent with these states being 
the J p = 0+ (D sJ (2317)) and J p = 1+ (£> sJ (2457)) 
members of the cs system (with L = 1, and "light 
quark angular momentum" j q = | ) , but their masses 
are a complete mystery. We thought we understood 
the cq mesons . . . but it appears that we don't. This 
topic will be covered further in Jussara de Miranda's 
talk on Standard Model charm studies. 5 



a In the case of the D sJ (2457), the Belle A, 7 results 4 ' 3 rule 
out J = 0, 2, and are consistent with J = 1. Note that an im- 
portant "exotic" hypothesis — that the D a j (2317) is wholly or 
partly a DK bound state, and the D 3 j (2457) likewise a D* K 
state — is also consistent with the 0+ and 1+ assignments. 



The largest discrepancy between theory and ex- 
periment in the charm sector is also off-limits, as no- 
one imagines that new physics is responsible. But 
we still don't understand why the measured cross- 
section 6,7 for e + e~ — > iprj c is an order of magni- 
tude larger than the NRQCD prediction; the inge- 
nious suggestion 8 of e + e~ — > 7*7* — ► ipip contam- 
ination has now been ruled out. 9 The ip cc fraction 
in e + e~ — ^ ip X is likewise far "too large", almost 
saturating ip production. We had thought that we 
understood cc production at this energy . . . but it's 
pretty clear that we don't. Tomasz Skwarnicki will 
have something to say about this, and other devel- 
opments in charmonia, in the next presentation. 

1.2. What I Am Talking About 

After discussing the problems of obtaining clean 
new physics signatures in the charm sector (Sec. 2), 
the largest part of the talk treats D° — D° mixing 
searches (Sec. 3); there are important new results 
on both y C p (Sec 3.3) and D° -> K + Tr~ (Sec. 3.4). 
The first CPT and Lorentz invariance violation study 
in charm has recently been published, and we dis- 
cuss it briefly (Sec. 4). The remainder of the talk is 
given over to searches for rare and forbidden decays 
(Sec. 5). There are new results from BaBar, Belle, 
CDF, CLEO, and FOCUS— including, as I say, a first 
observation — but let's not get ahead of ourselves. 

2. Finding Clear New Physics Signatures 

The difficulty with finding a clear signature of New 
Physics in the charm sector is this: it can be hard to 
know what the Standard Model (SM) prediction is. 
One way of thinking about the problem is to consider 
the masses of the quarks. 
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• The up and down quark masses are both small 
compared to the hadronic mass scale: m u < 
m d -C ^QCD- b Isospin is therefore a rather 
good symmetry, and (including now the strange 
quark) SU(3) of flavor, while broken, is useful. 

• At the other extreme, the beauty quark has 
m b > Aqcd, and can be considered as a high- 
energy physics "particle" : a "billiard ball with 
quantum numbers attached" . One can think in 
terms of the Feynman diagrams, in 6-sector pro- 
cesses, and not be too seriously misled. 

• The charm quark lies between the two extremes: 
m c > ^qcd, neither light nor truly heavy. 

The awkwardness of the charm mass thus puts lim- 
its on both symmetry- and quark-based thinking as 
guides to charm physics. If light hadron work is like 
swimming in the ocean, and 6-physics is like flying 
through the air, then in charm studies one is wading 
knee-deep through the brown muck. 

One should really think in terms of hadrons, not 
quarks, in charm. So-called "long-distance" con- 
tributions are important in many processes: quark 
loops are typically suppressed, so that hadronic pro- 
cesses take a leading role. These are usually difficult 
to calculate, especially as the charm mass lies in the 
resonance region. So if some parameter is supposed 
to be small, but observed to be large, one should be 
cautious before claiming new physics: perhaps the 
Standard Model contribution has just been miscal- 
culated. As discussed in Sec. 1.1, there have already 
been two major surprises in the last two years. 

3. D° - D° Mixing 

3.1. Mixing in the Standard Model and 
Beyond 

There are particular pitfalls in the interpretation of 
charm mixing searches. The SM box diagrams for 
mixing (e.g. Fig. 1) are doubly Cabibbo-suppressed 
and suffer from very efficient cancellations (the GIM 
mechanism): the expected mixing rate due to such 
processes is negligible. Since most new physics sce- 
narios introduce new particles that couple to the SM 
fields, they induce new loop diagrams such as Fig. 2, 

l We set aside the related fact that m u and are elements 
of the theory, rather than straightforward "observations" . 
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Figure 1. Box diagram for D° - D° mixing in the SM. 
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Figure 2. Box diagram for D° — D° mixing in a New Physics 
model with an extra down-type quark. 



with no (or a lesser degree of) cancellations. The 
result is an enhancement of the mass splitting of the 
D° — D° eigenstates, x = AM/T; hence the common 
statement that "D — D a mixing with measurable x 
is a signal of New Physics" . 

As discussed in the previous section, however, 
hadronic processes cannot be neglected. Final states 
common to D° and D°, such as KK, tttt, Kit and 
Kit, couple the two neutral Z?'s (Fig. 3); such con- 
tributions cancel in the SU (3)p limit, but to the 
extent that SU(3) of flavor is broken, they induce 
mixing. One might suppose that only the lifetime 
splitting parameter y = Ar/2r would be affected, 
as the intermediate states are real. The true situa- 
tion is more complicated. To the extent that quark- 
hadron duality holds, mixing can be estimated using 
the Operator Product Expansion: a recent study 11 
finds x ~ y ~ O(10~ 3 ). An alternative approach 12 
relying directly on hadronic intermediate states sug- 
gests that y may be as large as 0(1%). So as far as 
x and y are concerned, mixing provides a clean new 
physics signal only if x 3> y ~ 1CP 3 . 

k, % k, % 

d« I M D° 



K, 71, it, K 

Figure 3. Sample diagram for D° — D° mixing due to common 
hadronic final states. 
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Figure 4. Five parameters describing D° — D° mixing, and their relationship with quantities measured in experiment. 



3.2. Mixing Parameters & Measurements 

The full parameter space for mixing is more rich 
than (x,y): some important parameters are shown 
in Fig. 4, together with experimentally measurable 
quantities. The parameters are poorly-known: the 
strong phase difference Sk-k between D° — > K + ir~ 
and K~n + amplitudes is unconstrained, as is the dif- 
ference in particle content of the eigenstates |f 1,2) = 
p\D }±q\D ), A = (|p| 2 -M 2 )/(|p| 2 + M 2 ). But the 
news is not all bad. This year has seen the first 
measurements relevant to the CP violating phase 
</> = arg (qA(D° -> K-K+)/pA(D° K~K+)) 
(Sec. 3.3.6), as well as a major new analysis of 
D" -> K+n~ (Sec. 3.4). 

3.3. Mixing: Lifetime Difference, yep, 
and CP Violation 

The most popular measurement in recent years, how- 
ever, has been yep- Defined as the normalized life- 
time difference of the D°-D° CP eigenstates, it is 
typically measured using the non-eigenstate decay 
D a — > K~ir + as a convenient reference: 

= r(CP+) - r(cp-) _ t(d° -> K-TT+) _ 

VCP ~ T(CP+) + T{CP-) ^ r(D° -> K-K+) 

= y cos 4> + xA sin 



where the last relation holds for small values of the 
parameters. In the CP-conserving limit = and 
A = 0, so yep = y, as one would expect. In this limit 
yep is not a new-physics search parameter (since new 
particles are expected to affect x; Sec. 3.1) but a tool 
for measuring the level of mixing due to the SM. 

3.3.1. y C p: the FOCUS measurement (2000) 

Three years ago, the FOCUS collaboration mea- 
sured j/cp using a relatively clean sample of 10,000 
D° — > K~K + events, c and a K~ir + sample ten times 
that size. 13 Both inclusive and D*-tagged decays 
were used, under FOCUS-standard reconstruction, 
particle identification, and vertex detachment cuts; 
the result was obtained from a binned maximum- 
likelihood (ML) fit to the distributions of reduced 
proper time t' = (l—Nai)/ /3jc, where I, o~i are the D° 
decay length and its error, and N the minimum re- 
quired detachment of the production and decay ver- 
tices. 

The result was surprisingly large: j/cp = (3.42 ± 
1.39 ± 0.74)%, over 2tr away from zero. There was 
considerable excitement at the thought that charm 
mixing might be within our grasp, partly due to the 
usual association of mixing with new physics. But as 

inclusion of charge-conjugate modes is implied throughout, 
unless the context makes clear that they are treated separately. 
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Figure 5. BaBar yep analysis: 16 invariant mass distribution 
for data (points), projection of the fit (curve), and fitted back- 
ground component (shaded) for the four samples. 



discussed above, the most natural interpretation of 
percent-level yep would be a large lifetime difference 
parameter y, due to Standard Model effects. The 
new physics interpretation was possible, but some- 
what forced: mixing with i>j and large CP viola- 
tion, sin^ - 0(1), A - 0(1). 

3.3.2. y C p: Belle and CLEO (2002) 

New measurements have followed in short order, us- 
ing D° produced in e+e~ — > cc interactions at the 
B-factories. Belle 14 used a sample somewhat larger 
and cleaner than that of FOCUS, performing an 
unbinncd maximum-likelihood (UML) fit to inclu- 
sive D° — ► K~n + and K + K~ decay time distri- 
butions, and measuring yep = (— 0.5 ± l.O+o'.g) %• 
CLEO 15 used a smaller sample (the 9.0 fbT 1 CLEO 
II. V run), required a D*-tag, and added D Q — > 
7r + 7r~ to the usual modes; they measured yep — 
(-1.2 ±2.5 ±1.4)%. Both results are manifestly 
consistent with zero, and each other — and the fact 
that both are negative has led many to discount the 
FOCUS result. But it's worth noting that the aver- 
age yep from the three is positive, and ~ 1%. 



3.3.3. yep: BaBar (2003) 

This year BaBar has released a comprehensive yep 
measurement 16 based on 91 fbT 1 of data includ- 
ing both D*-tagged D° -» K~-n~ , K+K~, tt+tt" 
events, and inclusive D° — > K + K~ . As usual in 
e+e~ analyses, backgrounds are suppressed by a 
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Figure 6. BaBar j/cp analysis: 16 proper time distribution for 
data (points), UML fit projection (open histogram), and fit- 
ted background component (shaded) for the four samples. The 
points beneath each plot show the bin-by-bin differences be- 
tween the data and the fit, divided by the statistical error. 



center-of-mass momentum cut and vertex quality cut 
on the D°, particle identification (PID) cuts on the 
daughter tracks, and (for the D*-tagged samples) 
cuts on the D* -decay pion (the "slow pion"). The 
resulting D° samples are shown in Fig. 5. For each 
sample, the fitted mass distribution is used to deter- 
mine the event-by-event probability that a D° can- 
didate belongs to the signal, as opposed to the back- 
ground under the peak. This probability is then in- 
cluded in the likelihood function for each event in an 
UML fit to the proper-time distribution. 

These distributions, and the results of the fits for 
each sample, are shown in Fig. 6. In each fit, the as- 
sumed "underlying" distributions of both signal and 
background are convolved with resolution functions 
based on a sum of gaussians: most of these terms 
have widths of the form Sal, where a\ is the event- 
by-event proper time error, and S is a scaling factor, 
meant to account for deficiencies in modelling of the 
detector, etc.. This method is common to the earlier 
Belle 14 and CLEO 15 analyses and reflects a consen- 
sus on time-distribution fitting at the B-factories. 

A blind analysis was performed to obtain the 
mixing parameter: the weighted average over the 
four modes is yep = (0.8 ± 0.4_tg; 5 ) %, the most pre- 
cise measurement to date. 
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3.3.4. y CP : Belle, D* -tagged (2003) 

A new analysis from Belle, contributed to this 
symposium, 17 takes a different approach. The idea 
is to find a robust resolution function which does 
not rely on the estimated proper-time error: this al- 
lows binned ML fits to be used throughout, so that 
the goodness-of-fit can be explicitly checked. D*- 
tagged D° — > K~ir + and K + K~ events from the full 
Belle dataset of 158 fb -1 are used, subject to stan- 
dard reconstruction cuts, and the requirement that 
the proper time be well- measured. 

The D° lifetime, yep, and the parameters of the 
proper-time resolution function are all determined in 
a single simultaneous binned ML fit to the K~w + 
and K + K~ samples. The form of the resolution 
function is simple: a sum of five gaussians with a 
common mean, fixed relative normalizations, and 
floating widths. The gaussian widths for K + K~ are 
constrained to be the same as those for K~ir + , up 
to a single scale factor which is common to all terms. 
This parameterization has been studied using Monte 
Carlo (MC) data, and proves to be very stable: all 
values determined in a full decay-time fit match those 
fitted to the true resolution function, within their rel- 
ative errors (Table 1). 

Table 1. Belle j/cp analysis: 17 comparison of the parameters 
obtained in MC from a fit to the resolution function, using 
MC truth information (3rd column) and from the decay time 
fit, using reconstructed information only (4th column). The 
fractions of the five gaussian terms, fixed from the resolution 
function fit, are also shown (2nd column). 







fitted values (fs; except a) 


par. 


fraction (%) 


resolution fit 


lifetime fit 


o\ 


26.1 


95.1 ±1.3 


94.4 ±1.7 


o"2 


50.4 


177.0 ±2.2 


179.0 ± 1.2 


o"3 


19.8 


328.7 ±7.4 


328.2 ±2.2 


(74 


3.1 


675.7 ±24.9 


664.4 ±8.5 


CT5 


0.6 


2199 ±95 


2225 ±70 


^0 


[common shift] 


-1.51 ±0.22 


-0.95 ±0.54 


a 


[Kir -» KK scale] 


1.043 ±0.004 


1.042 ± 0.007 



The proper time distribution for data is shown 
in Fig. 7, together with the result of the binned ML 
fit; the confidence level is 94%. (All fits in the anal- 
ysis have an acceptable confidence level.) The fit 
returns a D° lifetime (from the K~it + sample) of 
412.6 ± 1.1 fs, consistent with the world average, 18 
and a mixing parameter yep = (1.15±0.69±0.38)%; 
the result is preliminary. 
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Figure 7. Belle yep analysis: 17 proper time distribution for 
data (histogram), binned ML fit (solid curve), and background 
component (dashed), for the K~tt + (upper) and K + K~ 
(lower) samples. 

3.3.5. t/cp-' summary of results 

For completeness, the j/cp results are summarized in 
Table 2. Note the dominance of the numbers from 
Belle (the tagged result, 17 as noted, is preliminary) 
and BaBar. The implications of these results were 
not addressed at this point in the talk, but in the 
Discussion at the end. We preserve this order here. 



Table 2. Summary of j/cp results. 



Experiment 


J/CP (%) 




E791 19 


0.8 


± 


2.9 


± 


1.0 


FOCUS 13 


3.4 


± 


1.4 


± 


0.7 


Belle, untagged 14 


-0.5 


± 


1.0 


± 


0.8 


CLEO 15 


-1.2 


± 


2.5 


± 


1.4 


BaBar 16 


0.8 


± 


0.4 


fO 

-0 


5 

4 


Belle, tagged 17 


1.2 


± 


0.7 


± 


0.4 



3.3.6. "yep + + "■■ D°-D° mixing and CPV 

In yep analyses based on D*-tagged samples, a CP- 
violation study comes "for free" : one can simply 
compare the lifetime of K + K~ events with different 
flavor D-tags, within the same analysis framework 
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used for j/cp- Belle 1 defines a parameter 

_ t(D -» ATA") - f -> M) 
r ~~ f (£> -> KK) + f -> 

w — A.y cos <ft — x sin 0, 

where the notation L for "effective lifetime" recog- 
nises that an exponential is being fitted to distri- 
butions which may not be strictly exponential. In 
the absence of CP violation in mixing, i.e. A = 0, 
the asymmetry parameter A r = —x sin <f), measuring 
the CP violating phase <f) = ^g^A/pA) due to the 
interference of decay and mixing. d The "AY"" pa- 
rameter of BaBar 16 is similar, differing by a factor of 
(f + yep)- Most systematic errors cancel due to the 
use of a common final state. The experiments find 
the following values, 

AF = (-0.8 ± 0.6 ± 0.2)% (BaBar) 
A T = (-0.2 ± 0.6 ± 0.3)% (Belle prelim.), 

consistent with zero; the measurements are statisti- 
cally dominated and will continue to improve for the 
life of the B-factories. 

Unlike mixing in general, CP violation associ- 
ated with mixing is a robust new physics signal: 20 
all charm mixing phenomena in the SM are domi- 
nated by the first two generations, so CP violation 
must be small. Even for x ~ y = 0(1%), we expect 
A r < 10~ 4 . So any significant non-zero measure- 
ment by this technique would be evidence of new 
physics contributing to mixing. One can imagine a 
scenario where both Standard Model and new physics 
processes lead to percent-level values of the mix- 
ing parameters (y and x respectively), and the new 
physics contribution leads to 30%-level CP violation: 
evidence for both SM mixing (via yep) and non-SM 
processes (via Ar) would emerge by the end of the 
B-factory era. 

3.4. Mixing: D° K+tv~ 

Anothcr new BaBar analysis, 21 of D° — ► K + w~ 
decays, has brought hadronic mixing analyses back 
into prominence. Pioneered by CLEO, 22 the sophis- 
ticated analysis method is sensitive to both mass- 
(x) and lifetime-splitting (y) of the neutral D eigen- 
states, and has been considered the technique of 

d The analysis is thus a close analogue of the B°/B° — > ipKg 
analysis in the 6-sector, measuring sin20i [= sin 2/3]. 



choice for e + e~ machines. The BaBar analysis has 
been gestating for some time — there were prelim- 
inary presentations (without final fit results) two 
years ago — and the related analysis at Belle is still 
underway, with only an intermediate result (the 
"wrong-sign rate" for Kir) in the public domain. 23 

3.4.1. D° — > K + n~ : the analysis method 

"Wrong-sign" hadronic decays such as D° — > K + tt~ 
occur via two paths: mixing D° — > D followed 
by Cabibbo-favoured decay D° — > K + n~ , and di- 
rectly by doubly-Cabibbo-suppressed (DCS) decay 
D° -> K+tt-. The DCS decay thus forms a back- 
ground to the mixing signal, and the two must be 
separated by reconstructing the decay time of the 
D° and exploiting the different time distributions: 
e _t for DCS decay and t 2 e~ l for mixing, where the 
proper time t is in units of the D° lifetime. The in- 
terference term between the DCS and mixing paths, 
which goes as te - *, cannot be neglected: 24 in fact 
it provides most of the mixing sensitivity, since the 
DCS rate is much larger than the mixing rate and 
the interference term is thus intermediate in size. A 
complication of the method is that this term is pro- 
portional, not to the lifetime difference parameter y, 
but to the quantity y' = y cos Sk-k — x sin Sxn which 
has been "rotated" by the strong phase difference 
Sk-k between the D° — > K + tt~ and K~n + decays. 

The subtle difference in time distributions (e~*, 
te~*, t 2 e~ l ) means that the time structure of back- 
ground events must be well-understood to avoid fak- 
ing a mixing signal. This is important as background 
levels are relatively high. The method used by 
CLEO, 22 which has been followed by both BaBar 21 
and Belle, 23 is to (1) tag the initial D flavor by re- 
constructing D* + — ► D°ir + ; (2) categorize the back- 
grounds according to their proper-time distribution; 
(3) measure their relative levels by fitting the data 
distribution in (M, Q) where M is the Kir mass and 
Q = M(K+ir-ir+) - M(K+tt-) - is the en- 
ergy release in D* + decay; 6 and (4) fix the back- 
ground levels in the fit to proper time. This is dif- 
ficult, but manageable, and CLEO reported limits 
at the 95% confidence level of \x n < 0.041% and 
-5.8% < y' < 1.0% based on this method. 



e BaBar uses 8m = M(K+7r _ 7r+) - M(X+7r"), which differs 
from Q by a constant. 
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Figure 8. BaBar D° -* K+n~ analysis: 21 (a,c) K+ it mass 
and (b,d) D*+ — D° mass-difference distributions, for (a,b) 
signal- and (c,d) background-dominated regions. Shown are 
the data (points), the projection of the fit (open histogram), 
and the fitted combinatorial (light), unassociatcd pion (dark), 
and doublc-misidcntification background (black). 
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Figure 10. BaBar D° — > K+tt~ analysis: 21 fit results and 
confidence intervals for the mixing parameters (x' 2 ,y'). 



3.4.2. D° -> K+n~: the BaBar analysis 

The BaBar analysis follows the CLEO model closely, 
but with a larger data sample (57.1 fb -1 ), lower 
backgrounds (due to BaBar's superior PID), and a 
data selection and fitting procedure finalized while 
remaining "blind" to the mixing results. The 
if + 7r~ distributions are shown in Fig. 8: background 
events are divided into combinatorial, truc-D°-plus- 
unassociated-pion, and double-misidentification cat- 
egories. This last type, where K~tt + is misidentified 
as tt~K + , is small but significant: see Fig. 8(b,d). 
Such events are retained to avoid any distortion of 
the other backgrounds due to targeted rejection cuts. 
The fit used to establish the background levels de- 
scribes the data well. 

The mixing parameters are then determined 
using unbinned, extended maximum-likelihood fits 
to the D° K+tt- ("wrong sign") and D° -> 
K~n + ("right sign") data. The likelihood terms 
for the signal and various background time distribu- 
tions are formed from underlying distributions (ex- 
ponentials or delta functions) convolved with event- 
dependent resolution functions similar to those of the 
yep analysis (Sec. 3.3.3); the event-by-event signal- 
and background-fractions are determined from the 
(rriKin 5m) fit. The results are shown in Fig. 9 for 
regions in (m^, Sm) dominated by signal (Fig. 9(a)) 
and background (Fig. 9(b)) events. 



Four overall fits allowing DCS decay only; DCS 
decay and CP violation (treating D° and D° sepa- 
rately); DCS decay and mixing, but no CP violation; 
and all three effects, are performed: the results are 
shown in Table 3. No evidence for mixing or CP 
violation is found. A slightly negative (and thus un- 
physical) value of x' 2 is preferred by the fit: this is 
taken into account in reporting the results. 



Table 3. BaBar D° — > K+n analysis: 21 parameters returned 
by the full fit to D°, D°, and combined samples. 



Fit case 


Para- 


Fit result (/10 


- 3 ) 




meter 


D° 


D° 


D° + D° 


Mixing 
allowed 


p(±) 
n ws 


3.9 
-0.79 


3.2 
-0.17 


3.6 
-0.32 




y<(±) 


17 


12 


13 


No mixing 


d(±) 

n ws 


3.9 


3.2 


3.6 



A rather careful procedure based on toy MC ex- 
periments is used to set frequentist confidence inter- 
vals in the fitting parameters x' 2 and y'\ the results 
are shown in Fig. 10. A remarkable feature of the 
analysis is that the allowed region in (a/ 2 , y') is com- 
parable in size to that of CLEO, despite the larger 
and cleaner dataset. Simulations by both BaBar 25 
and Belle show that for a given experiment, when 
the preferred value has y' > (as in the BaBar anal- 
ysis) , the allowed region becomes large compared to 
that when the data prefers y' < (as in CLEO's 
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Figure 9. BaBar D° — > K^ir~ analysis: 21 proper time distributions for (a) signal- and (b) background-dominated regions. Shown 
are the data (points), the projection of the fit (open histogram), and the fitted combinatorial (light), unassociated pion (dark), 
and double-misidcntification background (black). 



case) J It is thus almost meaningless to "average" 
the results of experiments whose preferred values fall 
in different y' regions. 

3.5. Mixing: Issues, Future Measurements 

A combined D° — > K + tt~ result will require a com- 
bined analysis: perhaps there is a need for a joint 
working group, or at least detailed consultation be- 
tween experiments, as we enter the CLEO-c era. 
Comparison of (x',y') and yep results is more com- 
plicated still. The strong phase difference 8kk, be- 
queathed us by the mischievous god of the color 
force, must first be measured: a significant shift 
might occur due to final-state interactions (FSI), 
shown to be significant in D — > hh decays by isospin 
analyses of CLEO 26 and FOCUS; 27 the latter study 
provides evidence for inelastic FSI. 

The default option is to wait for results from the 
CLEO-c run at the f/>(3770), as one of the analyses 
exploiting the coherent D°D° state promises an er- 
ror in cos Skit of order ±0.05. 28 ' 29 The other option 
is a complete measurement of the DCS D — > Kir de- 
cays, of which only D — > K + ir~ is currently known. 
CLEO has recently placed a limit on D+ -> K+ir ; 26 
a measurement is presumably within the reach of 

■^This may also be understood in qualitative terms using the 
proper time distribution for wrong-sign decays, 

e _t f Rd + VRoy't + — ^—t 2 j , 

where Rd is the DCS decay rate. For y' < there is a par- 
tial cancellation between the interference (oc te — ') and mix- 
ing (oc t 2 e~ J ) terms, and the fit becomes sensitive to small 
changes in both the x' 2 and y' parameters. For y' > 0, no 
such cancellation takes place. 



the other _B-factory experiments. Measurement of 
D + '° — ► K°ir + ' a rates relies on the measurement of 
D + '° — ► ii£7r + '° decays: the asymmetry with the 
corresponding K s mode is proportional to the inter- 
ference between decay amplitudes via K° and K°. A 
method for this measurement has been demonstrated 
by Belle, with a preliminary result for D° — > i^°7T . 30 

A promising new analysis method exploits the 
D° — > K s tt + it~ final state: Cabibbo-favoured (e.g. 
K*~ir + ) and doubly-suppressed (K* + n~) decays in- 
terfere, allowing measurement of the strong phase 
differences for the various resonant submodes; a 
study of time-dependence then yields the mixing pa- 
rameters x and y. This method has the advantages of 
superior scaling properties (the fit measures x rather 
than x 12 ) and sensitivity to the sign of the mass split- 
ting, in addition to the measurement of phases. It 
is, however, unproven. It has been championed by 
CLEO, who have measured the resonant substruc- 
ture of the decay; 31 a mixing study will presumably 
require the large samples available at the other B- 
factories. Those samples may also provide useful sen- 
sitivity from scmilcptonic decays, which have fallen 
out of favor, although there is an interesting unpub- 
lished L>° -> K^+£-p e analysis by CLEO. 32 

Fits for CP violating effects in mixing are now 
routine, and will increasingly become the main focus 
of study. By the next Lepton Photon meeting, the 
other major development will be CLEO-c analyses 
exploiting opposite-side tagging, geometric signal-to- 
background ratios, and a coherent initial state lead- 
ing to new mixing and CP violation observables. 
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4. CPT and Lorentz Invariance Violation 

Even more general analyses would allow for viola- 
tion of CPT symmetry, a manifest signal of new 
physics. While such studies have been carried out for 
kaons and B mesons, no CPT violation search had 
been performed in the charm sector until the recent 
FOCUS analysis. 33 Using D*-tagged D° -> R--K+ 
decays, they searched for indirect CPT violation 
parametrized by £ = (An — Aa2)/(Ai — A2), where 
A is the 2x2 effective Hamiltonian governing the 
time evolution of neutral D mesons, An are its diag- 
onal elements, and Xj its eigenvalues. The measured 
quantity is the time-dependent rate asymmetry 

CPTU T(D° -► K+tt-)+T{D° -» K-TT+Y 

which reduces to (Re(£)y — lm(£)x) Tt for xt, yt <C 
x and y are the mixing parameters discussed 
above. FOCUS finds (Re(^)j/ - Im(£)a:) = 0.0083 ± 
0.0065 ± 0.0041, with a 95% confidence interval 
-0.0068 < (Re(f )y - Im(£)a;) < 0.0234, consistent 
with zero. Their paper cites as an example the case 
where D° and D° mix with parameters (x,y) — 
(0.0,0.01): in this scenario, 0.68 < Re£ < 2.34. 

Within a formalism that allows for violation of 
Lorentz invariance, £ may depend on the vector mo- 
mentum of the studied particles, and on siderial time; 
the relation is a function of parameters 34 Aao y x,Y,z, 
where (X, Y, Z) is a non-rotating coordinate system. 
FOCUS fits for these quantities by measuring the 
CPT-violating parameter £ in bins of siderial time: 
we can summarize the (complicated) result by not- 
ing that the various \Aa^\ < O (l0~ 12 ) GeV at 95% 
confidence for the case (x, y, Sk-k) — (0.01, 0.01, 15°). 
This is to be compared with limits of order 10~ 21 
for the K°-K° system: the difference reflects both 
the size of the available data samples and the rela- 
tive strength of mixing in the two systems. The Aa M 
may in principle vary with flavor, so this measure- 
ment is important for completeness, even though the 
sensitivity does not compete with that for kaons. 
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Figure 11. CLEO D° -> 77 analysis: 36 D*+ - D° mass- 
difference distributions for (a) D° — > Tr°n and (b) D° — > 77. 

As a result, it can be difficult to calculate the SM 
expectation with certainty. The best new physics 
signals are those decays with extremely small pre- 
dicted rates, or special features which allow SM and 
non-SM contributions to be distinguished. 

The decay D — > 77, with an expected branch- 
ing fraction of order 10~ 8 , is thus a new physics sig- 
nal (Sec. 5.1); but D° -> (Sec. 5.2), whose 
predicted rate is five orders of magnitude smaller, is 
a more reliable one. 35 In the case of decays D — » 
M + £~ , the most robust new physics signal is not the 
decay rate, but the dilepton mass spectrum, which 
exhibits marked differences between SM and new- 
physics predictions (Sec. 5.3). And in each case, 
variant modes violating lepton flavor or number con- 
servation can be added "for free" to the analysis, in 
the spirit of a lamp-post search: there is no uncer- 
tainty on Standard Model predictions of zero. We 
consider each of these cases in turn, before treating 
decays D° —> Vj in the final section (Sec. 5.4). 



5. Rare and Forbidden Decays 

Moving to less exotic possibilities, flavor-changing 
neutral currents (FCNC) have not yet been observed 
in the charm sector. As with mixing, the SM parton- 
level loop contributions are subject to powerful can- 
cellations, and long-distance contributions dominate. 



5.1. D° -h> 77 (CLEO) 

Decays to photons have very small contributions 
from parton-level processes, but significant ones from 
the vector meson dominance (VMD) mechanism: 
B(D° -> 77) would be - 3 x 10" 11 if only short- 
distance mechanisms contributed, but is expected to 
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be of order 10~ 8 due to long-distance SM processes. 
This mode has not previously been studied. CLEO 
has recently used D* + — » D°t: + events in 13.8 lb -1 
of data to conduct a D° — > 77 search, normalizing to 
the D° — > 7r 7r° mode. 36 Under standard cuts, aug- 
mented with a 7T° — > 77 veto on the photons form- 
ing the D° — > 77 candidates, they accumulate fairly 
clean event samples and then fit the distribution of 
energy release Q from the D* decay, to keep differ- 
ences between the 7r°7r° and 77 modes to a minimum. 
The results are shown in Fig. 11: the agreement be- 
tween the data and the MC simulation is remarkable. 

A limit B(D° -> 77) < 2.9 x 1CT 5 is found at the 
90% confidence level, some three orders of magnitude 
above the SM prediction. There is thus some room 
for a new physics signal if future experiments can 
improve on CLEO's sensitivity. 

5.2. D° -> p+fj,- (CDF) 

The expected rate for D° — » is much smaller, 

O(10~ 13 ), whereas R-parity violating (RPV) Super- 
symmetry could lead to a branching fraction as high 
as 3.5 x 10~ 6 , just smaller than the previous ex- 
perimental bound. 35 The dimuon decay is thus a 
straightforward new-physics search mode. Using 
their upgraded detector and trigger system, which al- 
lows them to select a charm decay sample, CDF have 
conducted a search for this mode in the early Run II 
data. 37 With a fairly straightforward blind analysis, 
using D*-tagged events and D° — > 7t + it~ decay as 
a normalization mode, they observe no events 
and set a limit B(D° -> < 2.4 x 10" 6 at 90% 

confidence, improving on the previous bound by a 
factor of two. The 0(1) event background estimate 
relies on interpolation from the sidebands — events 
with true muon(s) dominate over the misidentifica- 
tion background — and will need to be better under- 
stood to significantly improve the limit. Presumably 
this is achievable, and improvements to this channel 
will depend on the progress of Run II data-taking. 

5.3. D -> hit (FOCUS) 

The only analysis of a "basket" of rare decay modes 
in recent times is by FOCUS, who have searched 
for decays D~}\ — > h^fifi, where ft. ± = 7r ± , i4T ± . 38 
For definiteness we will take D + — > 7r + /i + /i~ as an 
example. The predicted SM rate for this decay is 
O(10 -6 ), while (e.g.) the allowed R-parity violating 




0.5 1.0 1.5 

m„„ [GeV] 



Figure 12. Predicted 35 dimuon mass distributions M(/i+/x — ) 
for D+ — » 7r+/i+£t — . The solid line shows the sum of the 
short- and long-distance contributions in the SM; the dotted 
line, the R-parity violating contribution from SUSY at the 
level allowed prior to the FOCUS measurement — see the text. 

contribution, 35 15 x 10~ 6 , saturates the previous ex- 
perimental limit. 39 There would therefore seem to be 
potential for further restriction of RPV parameters, 
but not for observation of a new physics signal. 

The SM contribution, however, is dominated by 
the path D + — » tt + V — » 7r + /i + /i~, where the V are 
vector mesons. The predicted dimuon mass spec- 
trum (Fig. 12) thus shows pronounced peaks at the 
p and 4> masses, which dominate the SM rate. By 
contrast the spectrum for RPV is relatively flat, so 
that a new physics contribution comparable to or 
even below the SM contribution could be resolved 
by comparing the M(fj, + fi~) distribution of observed 

events with the various predictions. 

+ + - 
7C \l \l 

N 4 [ Sideband = 19 




1.7 1.8 1.9 2 2.1 

Mass GeV/c 2 

Figure 13. FOCUS D + — * 7r+/i+/x~ analysis: 38 the histogram 
shows selected events; the filled entry is in the signal region. 
The signal region (hatched) and the region used for sideband 
subtraction (cross-hatched) are also shown. 
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FOCUS is beginning to probe this region. A 
histogram of selected events for the D + — > tt + fi + n~ 
analysis is shown in Fig. 13, and results for all the 
modes are listed in Table 4. The analysis proceeds 
using standard FOCUS detached vertex, hadron- and 
muon-identification cuts; cut values are selected and 
background estimates calculated using a very careful 
"dual bootstrap" method, to minimise possible selec- 
tion biasses. No evidence is seen for any of the de- 
cay modes; existing limits are everywhere improved, 
in some cases by an order of magnitude. For the 
D + — > 7r + ^+^~ and Df — > tt + (i + [i~ modes, the 
sensitivity is approaching the SM prediction. 

Table 4. FOCUS D -* hU analysis: 38 measured limit on the 
branching fraction, SM prediction, 40 previous best limit, 39,41 
and expected CLEO-c sensitivity 29 for each mode. (Dg sensi- 
tivities are scaled from those of D+ and are not official CLEO- 
c numbers.) All entries are (/10~ 6 ). 



Mode 








FOCUS 


SM 


Prev. 


CLEO-c 


D+ -, 


K 


V 




9.2 


0.007 


44 


1.5 


D+ -» 


K 




V 


13 




120 


1.5 


D+ -, 




M~ 


M + 


8.8 


1.0 


15 


1.5 


D+ -» 






M + 


4.8 




17 


1.5 




K 


V 


> + 


36 


0.043 


140 


15 




K 




v 


13 




180 


15 




7T + 




M + 


26 


6.1 


140 


15 




7T 




M + 


29 




82 


15 



Further progress is expected at CLEO-c, whose 
sensitivities 29 are also shown: an improvement by a 
factor 3 ~ 8 is foreseen for the D + modes, reaching 
the SM expectation in the case of D + — > 7r + /i + /i~ 
and therefore restricting further the RPV contri- 
bution. Any signal in the lepton-number violating 
modes Df-, — » h~ , albeit unexpected, would of 
course be an observation of new physics. 



5.4. D° —> (fry, 4>tv°, cfirj (Belle) 

Finally we turn to radiative decays D° — > Vj, an- 
other vector meson dominance process in the Stan- 
dard Model. The Belle collaboration has conducted 
a search for D° — > (fry, exploiting double kaon iden- 
tification in (f> — > K + K~ to suppress backgrounds. 42 
Theoretical estimates for this mode, 43,44 dominated 
by D° — > UU' — > U7 (where the V" are vector 
mesons), are in the range (0.04 ~ 3.4) x 10 -5 , well 
below the previous limit of 1.9 x 10~ 4 but partially 
overlapping the sensitivity at the B-factories. 
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Figure f4. Belle D° — > (fry, <f>n°, <f>rj analysis: 42 <fiir° mass dis- 
tribution for data (points) and MC (histogram); the ^-mass 
sideband is also shown (shaded). 



D*-tagging and cuts on the D* momentum and 7 
energy are used to suppress the various combinatorial 
backgrounds. The dominant remaining background 
is due to the Cabibbo- and color-suppressed decays 
D° — > 07T° and <fir\, which have not previously been 
observed. With analagous cuts to select D° — ► (fin 
Belle sees a very clear signal in M(<pTr°) (Fig. 14) 
and the expected distribution of the hclicity angle of 
the <f> meson (not shown); the <f> is polarized in the 
D° decay. A smaller but still clear signal of 31 ± 9.8 
events is seen for D — > (fir], where a veto is imposed 
on photons for the rj — > 77 candidate consistent with 
belonging to a ir° — > 77 decay. 



-f data 
- best fit 
-- BG of best fit 
□ Feed-across BG 
D"->.<])7t n 

D°-»<j)r| 

D*-*(j) JI*7I° 




Figure f5. Belle D° — > c/yy, <frn°, <f>V analysis: 42 07 mass dis- 
tribution for data (points), the ML fit (open histogram), the 
background component of the fit (dashed), and the sum of 
D° —> <f>TT° , <f>r] and D + — > <f>ir + n° backgrounds (shaded). 



The contribution of these decays to the (fry spec- 
trum can then be reliably estimated; it is suppressed 
by a helicity angle cut |cos#hei| < 0.4, favoring the 
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Figure 16. Belle D° — > $7, 07r°, cAr; analysis: 42 cross-check 
of the <fyy result using the cos0 he i distribution. Left plot: 
data (points), and MC predictions for the total (open his- 
togram), total background (light), and non-07r° background 
(dark). Right plot: background-subtracted distribution for 
the data (points) and the MC prediction (histogram). 

transversely polarized of D° — > 07 over the lon- 
gitudinally polarized of the 07r° and iprj modes. 
The resulting <jyy invariant mass spectrum shows a 

clear D° -► 7 signal of 27.6±^;f jlj'-o events ( Fi S- 15 )' 
and is well-described by a linear combinatorial back- 
ground, the expected 07r° and 0?y contribution, plus 
the signal. The helicity angle distribution (Fig. 16) 
is likewise consistent with expectations. 
The final results of the analysis are 

B(D° -> 0tt°) = (8.01 ± 0.26 ± 0.46) x 10~ 4 
B(D° -► 077) = (1.48 ± 0.47 ± 0.09) x 10~ 4 
^ o -0 7 ) = (2.6O™°:^)xlO" 5 ; 

the 07 mode is the first radiative decay, and the 
first FCNC decay, observed in the D meson system. 
The measured branching fraction is at the upper end 
of the VMD predictions — consistent with Standard 
Model expectations — and is in no sense a new physics 
measurement, but it provides the first experimental 
reference point for predictions of other FCNC decays. 

6. Summary and Prospect 

While all observations are still consistent with Stan- 
dard Model expectations, there has been significant 
progress in charm analyses sensitive to new physics 
effects. The D-mixing measurement using D° — > 
K + n~ is now mature, with a major result released 
by BaBar, and a Belle analysis in the pipeline; the 
.B-factories will (presumably) exhaust this difficult 



technique. The sensitivity to yep continues to im- 
prove, and a 1% measurement is within reach of the 
current facilities during their projected life. Fitting 
for CP- violating effects has become standard, and as 
a robust test for new physics — and a technique still 
dominated by statistical errors — will become increas- 
ingly important in the field. And measurements of 
the new quantities observable at CLEO-c, and the 
promising but yet-untried D° — > Kgir + ir~ mixing 
analysis, are eagerly awaited. 

The first flavor-changing neutral current decay 
in the charm sector, D° — > 07, has finally been 
seen, and is consistent with Standard Model ex- 
pectations due to vector meson dominance. While 
not the most exciting possible channel — 77 or 
would have been respectively a shock, and a total 
revolution in the field — the 07 observation provides 
an experimentally-measured point where none previ- 
ously existed, and will presumably allow more precise 
SM predictions for other FCNC decays in the future. 
In the search for other rare and exotic processes, 
there are continuing improvements in the range of 
channels studied (with the first D° — ► 77 limit an- 
nounced, and the first test of CPT/Lorentz invari- 
ance conducted) and the reach of existing analy- 
ses (with the D -> Ml results from FOCUS). As 
with mixing, there will be significant contributions 
to these searches from CLEO-c in the next two years, 
and we all — including even the charm coordinators 
at BaBar, Belle, CDF, and FOCUS— are looking for- 
ward to new results in this new era. 
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DISCUSSION 

Hal Evans (Columbia): What is the average yep 
measurement? Or is there a reason not to cal- 
culate it? 

Bruce Yabsley: The average would be dominated 
by preliminary measurements from BaBar and 
Belle. Further, results appear to be "clumping" 
based on measurement technique indicating the 
possibility of systematic problems. That being 
said, the speaker's average is (yep) = (0-9 ± 
0.4)% 

Further Discussion, Added since LP2003: 

It's appropriate to return to this question in more 
detail for the written version of the talk. I reserved 
the yep average for the Discussion — assuming, cor- 
rectly, that someone would bring it up — partly for 
lack of time and partly for lack of a clear idea of how 
to treat it: any average is dominated by a prelim- 
inary number from Belle 17 and results from BaBar 
which, though "final" , had not been published at the 
time of the symposium. Since then, BaBar's paper 
has appeared in PRL, 16 so my reservations are now 
diminished. The situation, however, is still unclear. 



Table 5. Expanded summary of i/rjp results. 



Technique 


Experiment 


yep (%) 


Fixed 
target 


E791 19 
FOCUS 13 
My average: 


0.8 ±2.9 ±1.0 
3.4 ±1.4 ±0.7 
2.9 ± 1.4 




Belle 14 


-0.5 ±1.0 ±0.8 


e~^~e~ , 


CLEO 15 


-1.2 ±2.5 ±1.4 


untagged 


BaBar 16 


0.2±0.5i°;4 




My average: 


0.0 ±0.6 




BaBar, 16 K+ K~ 


1.5 ±0.8 ±0.5 




BaBar, 16 tt+tt" 


1.7 ±l.2±li 


.D*-tagged 


Belle, 17 K + K~ 


1.2 ±0.7 ±0.4 




My average: 


1.4 ±0.6 


Speaker's grand average: 


0.9 ±0.4 



The results from the various experiments arc 
shown in Table 5, sorted by the type of experi- 
ment and the method used: fixed target (E791 19 and 
FOCUS 13 ), e+e- with inclusive D° samples (Belle, 14 
CLEO, 15 and BaBar 16 ), and e+e" with L>*-tagged 
samples (BaBar 16 and Belle 17 ). The individual re- 



sults within the BaBar analysis have been listed sep- 
arately for this purpose. There is a clear clustering 
of the results according to technique: fixed target 
measures high ({yep) = 2.9%), e + e~ measures null 
{{yep) = 0.0%), and the D*-tagged analyses mea- 
sure an intermediate value ((yep) — 1.4%). These 
last results are completely dominant, as can be seen 
in Fig. 17, where the data are shown in the form of 
the "ideograms" used by the PDG. 18 





A all data 


fixed target / 


\ fixed target 


and e + e _ untagged / 





-8 -6 -4 -2 2 4 6 8 

y CP (%) 



Figure 17. PDG-style "ideogram" of the yep data in Table 5: 
each measurement (with mean xi and total error ±6xi) is 
represented by a gaussian with central value Xi, error 5x{, 
and area proportional to X/Sxi- The sum of all curves (solid), 
the sum of fixed-target and e+e - untagged data (dashed), and 
the fixed-target data alone (dot-dashed) are shown. 

Statistically speaking, the data are consistent, 
and the average is now 2a away from zero. But it 
seems to me a bit worrying that the different tech- 
niques don't actively corroborate each other. Put 
another way, having first become excited about the 
(positive) FOCUS result, and then rushed to dis- 
count it in the light of (negative) e + e~ results, I 
think we should be slow to interpret Z?*-tagged re- 
sults which split the difference. 

This question will become urgent after the next 
round of data-taking at the B-factories, if the central 
value stays at yep > 1% as the total error shrinks. 
One convincing cross-check would be to analyse D° 
decays to CP-odd final states such as Kg(p°,oj, 4>) 
and even K^rf'' at the £?-factories: for the same 
?/cp value, the shift in the lifetime has the opposite 
sign. I know of no such plans at this stage. 



